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We study the phase diagram of the ionic Hubbard model (IHM) at half-filling on a Bethe lat¬ 
tice of infinite connectivity using dynamical mean field theory (DMFT), with two impurity solvers, 
namely, iterated perturbation theory (IPT) and continuous time quantum Monte Carlo (CTQMC). 
The physics of the IHM is governed by the competition between the staggered ionic potential A 
and the on-site Hubbard U. We find that for a finite A and at zero temperature, long range anti¬ 
ferromagnetic (AFM) order sets in beyond a threshold U = Uaf via a first order phase transition. 
For U smaller than Uaf the system is a correlated band insulator. Both the methods show a 
clear evidence for a quantum transition to a half-metal phase just after the AFM order is turned 
on, followed by the formation of an AFM insulator on further increasing U. We show that the 
results obtained within both the methods have good qualitative and quantitative consistency in the 
intermediate to strong coupling regime at zero temperature as well as at finite temperature. On 
increasing the temperature, the AFM order is lost via a first order phase transition at a transi¬ 
tion temperature Taf{U, A) (or, equivalently, on decreasing U below Uaf{T, A)), within both the 
methods, for weak to intermediate values of U/t. In the strongly correlated regime, where the effec¬ 
tive low energy Hamiltonian is the Heisenberg model, IPT is unable to capture the thermal (Neel) 
transition from the AFM phase to the paramagnetic phase, but the CTQMC does. At a finite 
temperature T, DMFT-I-CTQMC shows a second phase transition (not seen within DMFT-I-IPT) 
on increasing U beyond Uaf- At Un > Uaf, when the Neel temperature Tjv for the effective 
Heisenberg model becomes lower than T, the AFM order is lost via a second order transition. For 
U A, Tn ^ t^/U{l — x^) where x = 2A/17 and thus Tn increases with increase in A/U. In 
the 3-dimensional parameter space of {Ult,Tjt and A/t), as T increases, the surface of first order 
transition at Uaf{T, A) and that of the second order transition at Un{T,A) approach each other, 
shrinking the range over which the AFM order is stable. There is a line of tricritical points that 
separates the surfaces of first and second order phase transitions. 

PACS numbers: 71.10.Fd, 71.30.-|-h, 71.27.-|-a, 71.10.Hf 


INTRODUCTION 

The Hubbard model is a paradigm for studying elec¬ 
tron correlation effects in metallic systems in condensed 
matter physics. It has played an important role in under¬ 
standing how electron electron interactions can give rise 
to many interesting phases, including insulating, mag¬ 
netic and super-conducting phases. It is worthwhile to 
explore whether there are interesting correlation effects 
in systems that would be band insulators in the absence 
of electron-electron interactions. Perhaps the simplest 
model in which one can carry out this exploration is an 
extension of the Hubbard model, known as the Ionic Hub¬ 
bard model (IHM), with a staggered on-site “ionic” po¬ 
tential A added in. In the recent past the IHM has been 
studied in various dimensions by a variety of numerical 
and analytical tools [IHIl- In one-dimension [H-Q it has 
been shown to have a spontaneously dimerized phase, 
in the intermediate coupling regime, which separates the 
weakly coupled band insulator from the strong coupling 
Mott insulator. In higher dimensions {d > I), this model 
has been studied by many groups using the dynamical 
mean field theory (DMFT) [I. IR 11^. IT^ . determi- 

nantal quantum Monte carlo [S Q: cluster DMFT Q 


and coherent potential approximation llj. The solu¬ 
tion of the DMFT self consistent equations in the para¬ 
magnetic (PM) sector at half filling at zero tempera¬ 
ture shows an intervening correlation induced metallic 
phase [1-0,0 [llj for intermediate values of the interac¬ 
tion U. When one allows for spontaneous spin symmetry 
breaking the transition from paramagnetic band insula¬ 
tor (PM BI) to AFM insulator generall y p reempts the 
formation of the para-metallic phase 0 [I^, except, as 
shown in a recent paper coauthored by two of us |I2l | us¬ 
ing DMFT with Iterated Perturbation Theory (IPT) as 
the impurity solver, for a sliver of a half-metallic AFM 
phase. Upon doping one gets a broad ferrimangetic half¬ 
metal phase [i^ sandwiched between a weakly correlated 
PM metal for small U and a strongly correlated metal for 
large U. 


In this paper, we provide a detailed discussion of the 
properties of the half-filled IHM on the Bethe lattice 
of infinite connectivity solved using DMFT, not only at 
T — 0, specially the half-metallic AFM phase, but also 
at the finite temperature at a level much more exten¬ 
sive than explored before [l^ . The DMFT equations are 
solved allowing the possibility of an anti-ferromagnetic 
(AFM) order. We show that at any finite T, the system 
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shows two phase transitions as the Hubbard U is tuned 
for a fixed value of the ionic potential A. As U increases, 
first the AFM order turns on via a first order phase tran¬ 
sition at Uaf followed up by a continuous transition at 
Un > Uaf at which the staggered magnetisation drops 
to zero continuously. As T increases, Uaf increases while 
Un decreases, due to enhanced thermal fluctuations, such 
that the range in U over which the long range AFM order 
survives shrinks. At a certain Tt^p, we have a tri-critical 
point that separates the lines of the first and the second 
order phase transitions (for fixed A). In the 3D parame¬ 
ter space oiU — A — T, there is a line of tricritical points 
separating surfaces of first and second order transitions. 
We believe that these features of the phase diagram of 
the IHM have not been pointed out earlier. 

Our results come from a detailed study of the IHM 
model using DMFT with two different impurity solvers, 
namely, the iterated perturbation theory (IPT) and the 
continuous time quantum Monte-carlo (CTQMC- imple¬ 
mented using TRIQS package [I^), which also allows us 
to explore what is the interaction regime where IPT, an 
approximate semi-analytic impurity solver, works well. 
We show that the zero temperature phase diagram ob¬ 
tained within DMFT-I-IPT matches well (see Fig. [^l, 
both qualitatively and semi-quantitatively, with that ob¬ 
tained from the DMFT-I-CTQMC (working at very low 
temperatures). In both the methods we find that at zero 
temperature, and for a finite A, the long range AFM 
order sets in beyond a threshold value of 17, which we 
denote Uaf, via a first order phase transition. For U 
smaller than Uaf the system is a correlated band insu¬ 
lator in which the gap in the single particle density of 
states (DOS) reduces as U is increased. Both Uaf and 
the jump in the staggered magnetization at the transi¬ 
tion point increase with A. Both the methods show clear 
evidence of the formation of a half-metal phase just after 
the AFM order sets in, followed by the formation of an 
AFM insulator (AFM I) on further increasing U. Note 
that the half-metal AFM phase is missed out completely 
in the Hartree-Fock theory. 

For weak to intermediate {U ^ 6t and thus moder¬ 
ately strong) couplings, where the interesting metallic 
and half-metallic phases are realized in this model within 
DMFT-I-IPT, there is a good quantitative consistency be¬ 
tween the CTQMC and IPT results. However, in the 
limit of extremely strong coupling, where energetically it 
is favorable to project out the double occupancies and the 
ionic Hubbard model maps onto the effective Heisenberg 
model at half filling, one starts seeing deviation between 
the CTQMC and IPT results. While DMFT-kCTQMC 
correctly captures the physics of the effective Heisenberg 
model (as was also shown in [13|), perhaps not surpris¬ 
ingly, IPT fails to do so. At any finite temperature, 
CTQMC shows two phase transitions as the Hubbard U 
is tuned. First, with increasing U, the long range AFM 
order sets in via a first order jump in the staggered mag¬ 


netisation nis at Uaf{T). On increasing U further, rus 
first increases, reaches a maximum and then starts de¬ 
creasing with [7, finally going to zero via a continuous 
transition at C/ = Un(T). This transition happens when 
T just crosses Tn where Tn is the Neel temperature of 
the effective Heisenberg model obtained at U = Un- 

Within DMFT-I-IPT, at any finite T, only the first 
phase transition is seen as the Hubbard U is tuned. 
Within IPT also the AFM order sets in with a first order 
jump at Uaf- However, as U increases further, the AFM 
order saturates to unity and never goes to zero again. 
Thus the second phase transition from the AFM insula¬ 
tor to paramagnetic phase at large U is not captured by 
IPT. 

Consistent with this, the thermal phase transition for 
the half filled IHM is correctly captured within IPT only 
for weak to intermediate U, but is correctly described 
by the CTQMC for strong correlations as well. For all 
values of 17 > Uaf, the thermal transitions to the para¬ 
magnetic phase shown by IPT are always first order. On 
the other hand, CTQMC shows a first order transition 
only for small values of U/t- For U ^ 2A, the staggered 
magnetisation falls to zero across a continuous transition 
with increase in T. From weak to moderately strong val¬ 
ues of U /t, the transition temperature increases with Ujt 
in both the methods. But for U > 2A, while the tran¬ 
sition temperature keeps increasing with U within IPT, 
it goes as UU/{U^ — 4A^) within CTQMC, following the 
physics of the Heisenberg model as it should. 

The rest of this paper is organized as follows. In Sec¬ 
tion I we present the details of the model and the cal- 
culational methods used. Section II describes in detail 
the T = 0 phase diagram of IHM at half filling within 
IPT and CTQMC. Here we see a good qualitative and 
quantitative consistency between the two methods for a 
large range of parameters. Section III has details of the 
finite T phase diagram within IPT and its comparison 
to that obtained using CTQMC. We focus specifically on 
the regime of extreme correlations, where CTQMC works 
well but IPT does not. We end this paper with conclu¬ 
sions and discussions. In appendices A and C we present 
detailed discussion on the nature of the phase transition 
involved while appendix B shows a comparison of results 
within the DMFT and the HF theory. 


MODEL AND METHODS 

The model we consider has tight-binding electrons on 
a bipartite lattice (sub-lattices A and B) described by 
the Hamiltonian 

H = -t ^ [c\^cja + h-c] +A^m-A^rii 

i£A,j€B,(7 i^A i£B 

~h U ^ ^ (1) 
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Here t is the nearest neighbor hopping, U the Hubbard 
repulsion and A a one-body staggered potential which 
doubles the unit cell. The chemical potential is chosen 
to be /r = C//2, so that the average occupancy per site is 
{{ha) + {ns)) /2 = 1 , corresponding to “half-filling”. 


Dynamical Mean Field Theory (DMFT) 


Here we study this model using the DMFT approach. 
The DMFT ap pro ximation is exact in the limit of large 
dimensionality [IJ, ll5j and has been demonstrated to be 
successful in understanding the metal-insulator transi¬ 
tion Bin in the usual Hubbard model, which is the 
A = 0 limit of Eq. ([T]). We focus in this paper on the 
anti-ferromagnetic sector of Eq. o, for which it is con¬ 
venient to introduce the matrix Green’s function 




Cscr (h, icUyj) 


-1 


( 2 ) 


where a, (3 are sub-lattice {A,B) indices, a is the spin 
index, k belongs to the first Brillouin Zone (BZ) of one 
sub-lattice, iuin = (2n -f l)7rr and T is the temperature. 
The kinetic energy is described by the dispersion Ck, and 
Cyi(B)a = A-l-^-S^(B)^(ia;n). Within the DMFT 

approach the self energy is approximated as purely lo¬ 
cal 0. Thus the diagonal self-energies Yiaa{i^^n) are 
k-independent and the off-diagonal self-energies vanish 
(since the latter would couple the A and B sub-lattices). 

The DMFT approach includes local quantum fluctu¬ 
ations by mapping 00 the lattice problem onto a 
single-site or “impurity” with local interaction U hy¬ 
bridizing with a self-consistently determined bath as fol¬ 
lows. (i) We start with a guess for mg and 

dn and compute the local Gaaii^n) = 
rewritten as 


/ OO 

dc 

-OO 


Po(e) 


Ca<t (i^n) 


( 3 ) 

where, for a = A{B), a = B{A), and po(e) is the bare 
DOS for the lattice considered (see below). We actually 
need to solve the problem for only one sub lattice and use 
the relations Gact(*w„) = —GBa{—i^n) and I]^o.(*w„) = 
U — Y,B(y{—i^n) valid at half filling, (ii) We next deter¬ 
mine the “host Green’s function” [IJ, [isj Q^aa from the 
Dyson equation Q^aai^^n) = G~l{iuin) + E«cr(*w„). (hi) 
We solve the impurity problem to obtain = 

Sqct [Qoarrii^n)] (iv) We iterate steps (i), (ii) and (hi) till 
a self-consistent solution is obtained. 


IPT as impurity solver 

We use as our “impurity solver” in step (iii) a general¬ 
ization of the iterated perturbation theory (IPT) 00 


scheme which has the merit of giving semi-analytical re¬ 
sults directly in the real frequency (w’*' = a; -I- i0+) do¬ 
main. The IPT ansatz S^(^^(a;+) = -i- Aaa'^a}{uj+) 

is constructed to be (a) exact for U/t 1, (b) exact 
tor t/U = 0, and (c) exact in the large w limit for 
all U/t, which imposes various exact sum rules. Here 
'^aa = Uuaa Is the HF Self energy with 



and 

3 poo 

E^^^(a;+) = [/^ / da [paa{ei)paa{-e2)paa{e3)] 

.^1 J-oo 

[/(ei)/(-e2)/(e3) + /(-ei)/(£2)/(-£3)] 

a;+ — ei -I- e 2 — ea 


This has the form of the second order self¬ 
energy with pacri^i) = -Im[^oa<T(e+)]/7r, where 

^Oaa(w+) = is the Hartree 

corrected host Green’s function and /(e) is the 
Fermi function. From condition (c) above we find 
that Aq^(j — ^- 0.^(1 / [n.QQ.( 7 (l riQo.a-)] with 

noacr = ~~ f- dui Im Goaa{‘jJ~^) ■ Note that at half 
filling, since hao- = 1 — tt-Bo-, Aaa- is same for both the 
sub lattices. For simplicity, here we present the results 
for the solution of the DMFT equations on a Bethe 
lattice of connectivity z ^ oo. The hopping amplitude 
is re-scaled as t —>■ t/^/z to get a non-trivial limit, and 
the bare DOS is then given by po{e) = \/4t^ — e^/(27r<^), 
which greatly simplifies the integral in Eq. 


CTQMC as impurity solver 

In this section we describe briefly the state-of-the-art 
impurity solver, the continuous time quantum Monte- 
Garlo (GTQMC) using the hybridisation expansion 
method 0 |, in the context of the IHM. The impurity 
model (IM) at site a corresponding to the IHM can be 
written as 

Him, a 'y ^ (Cfc ~ Sq,A)/^^^/j;qct + ^ ^ 

kcr k<7 

h.c] -b C/nat?^a 4 ,- (/-t-SaA)^4^Cao- (6) 

where Sa = 1(—1) for a = A(B). Him, a describes 
the “impurity” in sub-lattice a coupled to the bath 
of / fermions through the hybridisation term Vka- It 
is straightforward to show, within a Grassmann func¬ 
tional integral formalism, that we can integrate out the 
fermionic bath variables in the partition function for the 
IM. After this step the partition function at site a be¬ 
comes 



( 7 ) 
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where and coao- are Grassmann variables represent¬ 
ing the fermionic “impurity” degrees of freedom at a site 
belonging to the a sub-lattice, and Sa is the functional, 

/•/3 

= / drdr/cj^^ (T-)Soa<T (t “ T/)coaa M 

a "'O 

rP 

+ / drC7nat('r)nQ,4,(r) . (8) 

do 

Here (iwn), the host Green’s function at site a, is 
related to the hybridisation amplitude 14a via the rela¬ 
tion 


^0Q(T(^rCn) — 5a^ 4- /T ^acr(^^n) (9) 

where Aaa(*w„) = Efe ,,^„-,,[+Lg„(a)A the hybridisa¬ 
tion function. On the Bethe lattice of infinite connectiv¬ 
ity, the self-consistent hybridization function for the IHM 
is given by /S.aa{iwn) = t^GaaiiWn), giving a simple re¬ 
lation between host Green’s function and lattice Green’s 
function as 


cl (d) 


0— O 




Ct(u) 


T/? 


Otg, 


0 


• ■ — ip 

c/(d) h 


FIG. 1: A pictorial representation of a configuration gener¬ 
ated by the CT-HYB algorithm with one up-spin and one 
down-spin electron. The total length of the segment in the r 
space for which an electron with spin cr lives is L„, and 0-[\, 
is the total length of overlap (in r space) for which electrons 
with both 4 and spins are present. 


and 



(15) 


= iw„ + + fl- t'^Gaa{iOJn) (10) 

Hence Sa can be re-expressed as 


In our case the matrix Aa is block-diagonal in “up” and 
“down” spin labels. Then the above equation simplifies 
to 


r/3 

Sa = SZc + X] / drdT/c5„^(r)Aa^(r - r/)coaa(T/) 
a do 


where 


SZc ^ E - M + SaA^ 

rP 

C 0 atT('r) + ?7 / dr noa^{r)noai{T) (12) 
^0 


oo 


— = TT / 

Zoa ^ ^ k, 


a kc=0 


1 

/ drf...dr- / dr'^..dr-^ 

cr- Jo Jo 


detAaa{TTCOaAT-Z)claA'Cl)-COaa{TZ)claA'cZ))sJ^^) 


The CT-HYB algorithm generates “configurations” 
corresponding to the terms in Eq. (16) with weights 
proportional to their contributions to the partition 
function Z. One such configuration is shown in Fig. [T] 


The partition function, Za given by Eq. [7] can then be 
expanded as a power series in S'S, 


^a = -^0aE^ / dri...dTfe / dT'i...dT4(ie<Aa X 
^ 0 0 

{TrCOaa, (Tl)cEq (r'l) ■ ■-CQaa^ (T-fe)cEq (^'fe)) cc (13) 


where 


Aa 


('^1 5 l) 

• Aaaiaa'k„i'Cl,T k)'\ 


■ Aaakaa'ki'Cki'T k) J 


( 14 ) 


The CT-HYB algorithm can calculate important quan¬ 
tities such as the finite temperature imaginary-time 
Greens function, the density, the double occupancy etc. 
For example, the occupancy Uaa is estimated from the 
average length of all the segments: naa={Laa)MC //3; the 
double occupancy is obtained from the overlap Ogti of 
segments as Da = {Oa\])MC/P', etc. For details see |l8l| . 

The DMFT self consistency loop run as follows. (1). 
One starts with a guess for the local Green’s function 
Gaaiiwn) where for a = A,B, a = B,A. (2). The host 
Green’s function for the a sub-lattice, Qoaaiiwn), is cal¬ 
culated using Eq. (fTOl) . (3). Using the host Green’s func¬ 
tion Goa-ai'iwn) the impurity solver calculates Gaaiiwn)- 
Then step 2 is invoked again, and the process is repeated 
until GA(B),cr{i'Wn) converges. We implement CT-HYB 
using TRIQS package [l9| . 
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FIG. 2: Phase diagram of the model in Eq. o at half filling 
obtained using DMFT for the Bethe lattice with IPT at T = 0 
(solid lines) and CTQMC done at T = 0.02t (dashed lines). 
A first order transition takes place at Uap such that for U > 
Uaf the system has long range AFM order, while for U < 
Uaf it is a PM BI. For Uhm > Uaf, the spectral gap in 
one of the spin components vanishes, resulting in a HM AFM 
phase for U = Uhm- For larger values of U the system is 
an AFM insulator. Note that the transition points obtained 
using the two methods are in fairly good agreement with each 
other. 

T = 0 PHASE DIAGRAM OF THE HALF-FILLED 

IHM 

The zero temperature phase diagram of the half-filled 
IHM obtained from the DMFT-I-IPT study and the 
DMFT-bCTQMC study (at T = 0.02t) is shown in Fig. [51 
With increasing U there occurs a first order transition be¬ 
tween the PM BI and an AFM phase, characterized by a 
non-zero staggered magnetization TOs , at some threshold 
U = Uaf (which is an increasing function of A). Inside 
the AFM phase, a half metal (HM) phase appears at 
U = Uhm > Uaf, where the gap in the single particle 
density of states (DOS) vanishes for one spin component 
while the other spin component has a non-zero spectral 
gap. When U increases well above Uhm, the system be¬ 
comes an AFM insulator (AFM I), where the gap in the 
DOS for both the spin components is controlled by, and 
increases linearly with, U. 

The phase diagram in Fig.|5]has been obtained from an 
analysis of various physical quantities, which we describe 
in detail below. 

Single particle density of states (DOS): In this 
subsection we discuss the single particle DOS /9a,cr(w) = 
~ Sfc Gaa{k, 0 J^)lTT, Calculated using DMFT-I-IPT. 
Here a represents the sub lattice A, B and a is the spin. 
Since at half filling pAa-i^) = PBa{—<^), we will focus 
only on the total DOS Pct(w) = pA(j{<jj) + Fig- El 

shows how P(t(w) evolves as a function of U for a fixed 
A = l.Ot. At small U < Uaf (=3.0t for A = l.Ot), there 
is spin-symmetry in the DOS, and Pcr(w) has a finite gap 



-1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 


P(M) P(w) 

FIG. 3: The single particle DOS, pc,{uj), calculated within 
DMFT-I-IPT, plotted as a function of uj for different values 
of U/t for A = l.Ot at n = 1. The red curves are for the up 
spin component and the blue dotted curves are for the down 
spin component. For U < Uaf = 3.0t, the DOS is same 
for both the spin components, with a non zero spectral gap 
which decreases as Ujt increases, and the system is a PM BI. 
For U > Uaf, the DOS becomes different for the two spin 
components. At U = 3.It > Uaf, the DOS for the up-spin 
component has a vanishingly small gap while the down spin 
DOS still has a finite gap. This is in close vicinity of the HM 
AFM point Uhm- At larger U values, there is a finite gap in 
the DOS for both the spin components, and the system is an 
AFM I. 

which decreases as Ujt increases, as shown in greater 
detail and clarity in Fig. [H We call this phase a PM 
BI as it is adiabatically connected to the U — 0 band 
insulator. 

For U > Uaf, the spin symmetry in the DOS is lost as 
seen in the top two plots of Fig. [31 The spectral gap in 
the up-spin component of the DOS is smaller than that 
for the down spin component, as can be seen more clearly 
in the inset of Fig. |4l We note that at half filling, even 
in the symmetry broken phase there is no net moment, 
i.e., n-\- = nj, = 1/2. This is because of the symmetry 
relations of the Green’s function (discussed earlier) which 
implies that nAa- = ^ and thus the total density 

of particles with spin cr is Uo- = ^[n-Aa + bbo] = 1/2, 
although from the top two plots of Fig. [3] it might seem 
that there is a net moment. The point is that Fig. |3| 
shows only the low ui DOS, where the area under the 
DOS for the up-spin component is larger than that for 
the down spin component due to smaller spectral gap for 
the up-spin component. However, the weight loss for the 
down spin component in the low uj regime is compensated 
by its large w part and the condition for no-net moment 
Uct = 1/2 holds. For U > Uaf, what the system has is 
a staggered moment, m-s = riAt — bai = bb\. — bb^, as 
discussed in more detail below. 

As U increases above Uaf, the gap in the up-spin com¬ 
ponent of the DOS decreases rapidly, and becomes van¬ 
ishingly small at a critical value U = Uhm (equal to 
3.09t when A = l.Ot), while the down-spin component 
still has a finite spectral gap (See Fig. 2]). Thus the IHM 
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FIG. 4: The spectral gaps Egf and Eg_i, obtained from the 
DOS within DMFT+IPT, plotted as functions of U for A = 
l.Ot at n = 1. The red points are for the np-spin component 
and the blue points are for the down-spin component. For 
U < Uaf in the BI phase, Eg^ = Eg^ and both decrease with 
increasing Ujt. At [/ = Uaf, there occurs a jump separating 
the two gaps, snch that Eg-^ is less than Eg^,. Eg^ becomes 
vanishingly small (< O.Olt) at U = 3.It, close to the HM 
AFM point. Both Egf and Egi increases with increase in Ujt 
in the AFM I phase {U > Uhm)- The inset shows Eg in the 
vicinity of the transition point. 


has a half-metal phase at a quantum critical point Uhm 
embedded within the AFM regime. This prediction is 
further reinforced below from the low w analysis of the 
spectral function. As U increases further, the spectral 
gap in the DOS opens up again for the up-spin compo¬ 
nent as well, with both the spectral gaps increasing with 
Ujt. This is the AFM insulating (AFM I) phase. 

Staggered magnetization and staggered occu¬ 
pancy: The staggered magnetization mg, defined as 
rus = (ttIzA — tnzB)l^, calculated both within the 
DMFT-bIPT (T = 0) and DMFT-bCTQMC (at T = 
0.02t) is shown in Fig.[5l For a given value of A, the stag¬ 
gered magnetization is zero below the corresponding 
U af and becomes nonzero for larger U , with a discontin¬ 
uous jump at Uaf corresponding to a first order phase 
transition between the PM BI and the AFM phase. Note 
that in the presence of the staggered potential, which 
opens up the gap in the DOS characteristic of the BI 
phase, the AFM instability does not occur unless U ex¬ 
ceeds a finite thresh hold value Uaf- The larger the value 
of A, the larger is the value of U required to overcome 
the effect of A and turn on the magnetisation. Thus both 
Uaf and the jump in at Uaf are increasing functions 
of A. The bottom panel of Fig. [5] shows the staggered 
occupancy, i.e., the difference in filling factor on the two 
sub lattices, defined as Sn = {ub — ua)!"^- Due to the 
staggered on site potential, this difference is always non 
zero, even though the Hubbard U tries to suppress it. 
For U < Uaf, Sn decreases monotonically and rapidly as 
a function of U. At Uaf, there occurs a discontinuity in 
Sn. For U > Uaf, Sn decreases more slowly with increas¬ 



FIG. 5: Top panel: Staggered magnetization rUs plotted as a 
function of Ujt at half-filling. A first order phase transition 
takes place with the onset of ma at Uaf- Bottom Panel: 
Staggered occupancy 5n plotted as a function of Ujt at half¬ 
filling. dn is non zero for all values oiUjt and a discontinuity 
occurs in in at Uaf- In both the panels the points connected 
with solid lines represent the data obtained from DMFT-I-IPT 
at T = 0 and the points connected with dashed lines show 
data obtained within DMFT-I-CTQMC at T = 0.02t. There 
is quantitative consistency between the two methods for a 
range of A values. The phase transition is clearly first order 
in both the methods. 


ing U, but eventually becomes rather small in the AFM 
I phase, asymptotically approaching zero (as tA/C/^ ) as 
U —>■ oo. 

Note that for all the A values, nig obtained using 
the CTQMC solver is slightly smaller than that from 
the IPT solver, while the transition point Uaf obtained 
using CTQMC is larger than that within IPT. This is 
because CTQMC captures the effects of quantum fluc¬ 
tuations better than IPT. But overall, in the small to 
intermediate Ujt regime, there is good quantitative cor¬ 
respondence between the low temperature CTQMC data 
and the T — 0 data obtained within DMFT-I-IPT. Also, 
the nature of the phase transition is the same in both the 
methods. 

The results in Fig. [5] have been obtained by solving 
the DMFT-I-IPT equations starting from a small U value 
and increasing U slowly. When the DMFT equations are 
solved starting from a large U guess and then decreasing 
U slowly, one gets a different curve for mg (and also for 
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(5n) (Fig. [131). A comparison of the ground state energies 
of these two spin-asymmetric solutions for the DMFT 
equation with the ground state energy of the PM sector 
shows that the real transition point Uaf is the one where 
rus becomes non zero for the first time coming from the 
small U side. The hysteresis analysis discussed in Ap¬ 
pendix A for A = l.Ot confirms the nature of the transi¬ 
tion from the PM to the AFM phase as being first order. 
But for very small values of A, where both the transition 
point Uaf and the jump in magnetisation at the transi¬ 
tion point are very small, numerically it is difficult to see 
the nature of the transition. Since in the small U regime, 
the Hartree-Fock (HF) theory also works well (as shown 
in Appendix B), we have carried out Ginzburg-Landau 
(GL) expansion of the ground state energy within the HF 
theory and confirmed that the phase transition from the 
PM to the AFM phase is of first order for any non zero 
A (for details see Appendix G). 



2.0 2.5 3.0 3.5 4.0 4.5 


U/t 

FIG. 6: F-i- and Fj, (see Eq. 18) plotted as functions of 
U/t for A = l.Ot. Points connected with solid lines are ob¬ 
tained within DMFT-I-IPT (T = 0) and points connected 
with dashed lines are obtained within DMFT-tCTQMC (T = 
0.02t). Note that F^ changes sign within the AFM phase for 
U > Uaf and crosses zero at Uhm = 3.09t within IPT and 
Ufm = 3.25t within CTQMC for T = 0.02t. 

Low uj analysis of the spectral function: To un¬ 
derstand the trend of the spectral gap and to confirm 
the existence of the HM AFM phase, we have carried out 
a low w analysis of the self energy and the single parti¬ 
cle spectral function. The IPT self energy = 

S'^^(cl’) + zS'/^(cj) has S"(a;) vanishing for |w| < 3Eg^ in 
both the insulating phases. This can be understood from 
the imaginary part of of Eq. ®, which comes from 
a three fermion final state. However, this is an artifact 
of ignoring collective modes (spin waves) within DMFT. 
In reality, since there are gapless spin waves that can 
be excited, the imaginary part of self energy will be zero 
only for |a;| < Egu and the phase space constraints would 
make the result for just above threshold quite small. 
In the discussion below, we assume that E"^(w) = 0 for 
||^| ^ Egtr- 


In both the insulating phases, can be writ¬ 

ten at low w as a Taylor expansion E(j^(a;) = E(j,g.(0) -I- 
w -|- ..., where can be shown to be in¬ 
dependent of a. The spectral function is defined by 
Aaaa{e,oj) = (-I/tt) ImGaca(c,w+). Since E"^ = 0 
for |a;| < 3Ega, we find from Eq. that Aaaai^,^) = 
S{rcr{u}) — e^) with r^{u}) = {uj + fi — A — E'^^{u})){uj + fj.+ 
A — E^^(a;)). As e is real, w’s which satisfy ra{oj) < 0 
lie within the gap. The energy gap is then given by 
ra- (Ego-) = 0 which, using the low-energy form of 
given above leads to the result 

Eg, = Z,\A-U/2 + E'^,{uj = 0)\ 

= Z,\A + Ul2-E'j,,{uj = 0)\ (17) 

where we have used the particle-hole symmetry. Let 
us write E'^^{ui = 0) = Saa + Uuaa where the second 
term on the right hand side is the self energy within the 
Hartree-Fock approximation. Then one gets a more elab¬ 
orate form for the expression of the gap, which is given 
below: 

Ega = Za\A - U/2{5n + anis) + 

= Z,\A - U/2{5n + ams) - SbA = Z,\T„\ (18) 

Fig. |6| shows Fcr as a function of [/ for A = l.Ot ob¬ 
tained within DMFT-tIPT and DMFT-tGTQMC(r = 
0.02t). Within the CTQMG, the Green’s function can 
be calculated only at Matsubara frequencies, and thus 
the single particle DOS and the spectral gaps can not be 
obtained directly from the CTQMG data. But the low 
energy part of the self energy can be used to get an es¬ 
timate of the spectral gap even from the CTQMC data. 
Specifically, in the CTQMC data, we have extrapolated 
the self energy to zero frequency and obtained approxi¬ 
mate values of Saa- 

For U < Uaf, in the PM BI phase, F^ = F^ > 0 
and decreases as U/t increases for a given A. At Uaf, 
F^ becomes different from F^. As U/t increases further, 
within the AFM phase, F^ decreases and becomes neg¬ 
ative for U > 3.It within IPT and U > 3.25t within 
CTQMC. Thus it must pass through a zero, making Eg-f 
zero inside the AFM phase, at Uhm = 3.09t(3.25t) for 
A = l.Ot within IPT(CTQMC). On the other hand, 
F^ remains always positive, giving a non-zero spectral 
gap for the down-spin component for all values of U/t 
including Uhm- Note that Za is always positive and 
less than one by definition, and that F^ obtained from 
CTQMC and IPT show good quantitative correspon¬ 
dence for U < Uaf- As U increases further, iFo-l within 
IPT becomes much larger than that within CTQMC. 
Within both the methods, we do see a novel, half metal 
AFM phase at Uhm > Uaf, inside the AFM phase of 
the correlated BI. Note that the half-metal AFM phase is 
missed out completely in a simple mean field theory like 
Hartree-Fock theory, though the BI to AFM transition is 
captured (see Appendix B). 









FIG. 7: The kinetic energy {IC^} and (/C 4 .) plotted as func¬ 
tions of U/t for A = l.Ot. The points connected with solid 
line are obtained within DMFT-hlPT (T = 0) and those con¬ 
nected with dashed line are obtained within DMFT-I-CTQMC 
(T = 0.02t). The kinetic energy decreases with increase in 
U/t deep in the PM BI phase, which indicates states becom¬ 
ing more extended with increasing U, while it increases in the 
AFM I phase, suggesting localization. For Uaf < U < Uhm, 
(/C-f) decreases with U/t, reaching a minimum at Uhm, while 
(Ai 4 ,) increases with increase in U/t as in the AFM I phase. 


FIG. 8 : Average double occupancy Da = {na^Uai) on sub¬ 
lattice a = A,B vs U for A = l.Ot. Points connected 
with solid line are calculated within DMFT+IPT(r = 0) 
and points connected by dashed line are calculated within 
DMFT-I-CTQMC (T = 0.02t). Due to the staggered po¬ 
tential, Da ^ Db for all values of U/t with Db showing 
a monotonic decrease with U/t. Note that Da within both 
the methods matches fairly well for U < 6t. 


Kinetic energy; We have also studied the spin- 
resolved kinetic energy which is defined as (/Co-) = 
“I I-oo^^ / e Po(e)-fwG'^ 5 (e,w+). Fig. [7] shows 
the (JCa) obtained within DMFT-fIPT (T = 0) and 
DMFT-hCTQMC (T = 0.02t). To calculate the 
KE within DMFT-I-CTQMC, which gives the Green’s 
function at fermionic Matsubara frequencies, we use 
TY.n^AB{i^n) = - ^ j duj f {u:)XmG’/^g{u)+), where 
/(w) is the Fermi distribution function, and derive the 
following version of the above expression for the KE: 

(/Co) = 2r f dee po{e) ^ G^B(e, /a;„) (19) 

n 

where G^^(e,ia;„) is the off-diagonal element of the full 
Green’s function defined in Eq. |31 

In the PM BI phase, as the spectral gap reduces with 
increase in U, (/Co) decreases until the correlation starts 
pushing the spectral weight from low energy to higher 
energy region. Once this happens, even though the spec¬ 
tral gap is decreasing within the BI phase, there occurs 
a slight increase in (/Co). 

In the AFM I phase, the kinetic energy for both the 
spin components increases with increase in U due to the 
increase in the spectral gap. In the regime for Uaf < 
U < Uhm, (/ Ct ) decreases with increase in U just like 
in small U limit of the BI phase. On the other hand, 
(/Cj,) starts increasing with U like in the AFM I phase. 
Note that the kinetic energy for the up-spin component 
is minimum at Uhm where the spectral gap is zero for 
the up-spin component and we have a HM AFM. 


Double occupancy: The average double occupancy 
at site a, Da = {na^na]) can be calculated using the 
following equation within IPT: 


Da — {Ba'\B.a\/} 


1 


[TY^i^nGaA 

71,(7 


iuJn) + Pan 


a 


(ic)] 


( 20 ) 

with a = B,A and pa = (z-^ + SaA). Within CTQMC we 
calculated Da by directly calculating the trace of na^Uai- 
Fig. m shows Da for A = 1.0/ at T = 0 obtained within 
IPT (T = 0) and CTQMC at T = 0.02/. For the IBM, 
since a non zero A// prefers to put more holes on the A 
sublattice and more double occupancies on the B sublat¬ 
tice, for all values of U/t, Da 'C Db- As U/t increases 
Db shows a monotonic decrease with a discontinuity at 
Uaf- Da on the other hand, first increases slightly as U 
increases below Uaf and then starts decreasing with U. 
Fig. |8] clearly shows that up to moderately strong values 
of U/t, the average double occupancy within IPT is quan¬ 
titatively very close to that obtained within CTQMC. 

All of the above analysis shows clearly that for the 
U/t range from weak to moderately strong, the IPT and 
CTQMC results match well. In our discussions in the 
following sections we focus on the differences between two 
approaches that arise when one looks at the extremely 
correlated regime of the IHM. 


EXTREMELY CORRELATED REGIME OF THE 
IHM AND FINITE T PHASE DIAGRAM 

In this section we consider the extremely correlated 
regime of the IHM, namely, U ^ t, A. In a regular 
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Hubbard model (A = 0), the limit of 17 ^ t effectively 
projects out doubly occupied sites from the Hilbert space. 
For the IHM, at half filling, the energy cost for having a 
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FIG. 9: Top: The staggered magnetisation vs Ujt for T = 
2t/35 and A = l.Ot calculated within DMFT+CTQMC. As 
U increases, first the AFM order turns on at Uaf with a 
jump in rris. As U increases further, at Un Uaf, the 
AFM order goes to zero continuously. Bottom: The staggered 
magnetisation rus vs U/t for T = 0.4t and A = l.Ot calculated 
within DMFT+IPT. As U increases, tUs keeps increasing and 
finally saturates to its maximum value. 


DMFT+IPT 



double occupancy on A(B) sublattice is 17 ± 2A. Thus, 
it is only for 17 t and U 2A that one obtains the 
extremely correlated regime of the IHM where doubly 
occupied sites are projected out. In this limit, at half¬ 
filling, the effective low energy Hamiltonian for the IHM 
is again the Heisenberg model 

^eff = • Sj — ninj/4:], ( 21 ) 

(b> 


but the spin-exchange coupling is now J = J/{1 — x^) 
with X = 2A/17 and J = /U. The Neel temperature 

Tjv of the Heisenberg model is proportional to J, which 
therefore depends upon 17. In dimensions higher than 2, 
for temperatures lower than Tm, the system has AFM 
ordering but the order is lost via a continuous transi¬ 
tion as T increases past Tjy. In a finite T calculation 


for the half-filled IHM, as U increases beyond 2A, J re¬ 
duces and eventually at some Un where the correspond¬ 
ing Tn{U = Un) > T, the magnetization is lost. Thus 
at any finite T, as one increases 17, two phase transi¬ 
tions should be seen for the half-filled IHM. First at Uaf, 
where the magnetisation sets in via a first order transi¬ 
tion, typically, and then at a larger Un > Uaf where the 
magnetisation is lost via a continuous phase transition. 


DMFT+IPT 



T/t 



FIG. 10: Top: The staggered magnetization rria vs T/t 
for A = l.Ot and various values of U obtained within 
DMFT-I-IPT. Bottom: ms vs T/t obtained within GTQMG 
for A = l.Ot for various values of U. Note that within IPT rris 
drops to zero via a first order phase transition at Taf which 
increases monotonically with U/t. But within GTQMG, for 
small values of U, though ms goes to zero via a first order 
transition but for larger values of U/t there is a clear continu¬ 
ous transition as a function of T in contrast to the IPT results. 
Also the transition temperature within GTQMG shows an in¬ 
crease with U only upto U = 5.Of and starts decreasing with 
further increase in U/t. 


This is exactly what we see in the CTQMC result as 
seen in Fig. [HI which shows the staggered magnetization 
vs 17/1 for A = 1. Of and T = 21/35. We see that the stag¬ 
gered magnetisation turns on with a first order phase 
transition at Uaf, increases with increasing U initially 
and then starts decreasing with further increase in 17, fi¬ 
nally vanishing a,t Un Uaf- But this second phase 
transition is not captured within DMFT-hlPT, for which 
the finite T the phase diagram is basically similar to the 
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T = 0 phase diagram. Once the AFM order sets in at 
Uaf, as we keep increasing U beyond Uaf, the staggered 
magnetization keeps increasing and never becomes zero, 
as shown in Fig. 13 Thus though the suppression of dou¬ 
ble occupancy for large U is captured correctly to some 
extent within IPT, spin physics and the physics of the vir¬ 
tual hopping resulting in the effective Heisenberg model 
is not captured. Hence, while IPT interpolates between 
weak coupling to strong coupling regime (by satisfying 
the atomic limit), at the end it is a second order per¬ 
turbation theory and especially for issues that crucially 
involve spin physics its validity breaks down in the regime 
of extremely strong correlations. 

THERMAL PHASE DIAGRAM 

Finally, we discuss how the AFM order is lost as the 
temperature T/t increases for a fixed value of A and U. 
The top panel of Fig. (TU] shows the finite temperature 
results obtained within DMFT-fIPT for A = l.Ot and a 
few values of U/t. As shown here, the staggered magneti¬ 
sation goes to zero via a clear first order phase transition 
at Taf- On the other hand, as shown in the bottom 
panel of Fig. [TOl within CTQMC the AFM order goes to 
zero via a first order transition, as the temperature T/t 
increases, only for small values of U/t. For U ^ 2A, 
the AFM order is lost continuously with a second order 
phase transition at Tjv. 

Further the transition temperature from the AFM 
phase to the PM phase has a very different dependence on 
U and A within IPT and CTQMC, specially for U ^ 2A. 
Within IPT, the transition temperature increases with 
increase in U for a fixed A irrespective of whether we are 
in the intermediate coupling regime or in the regime of 
extreme correlations. To be more specific, it follows U, 
and does not follow J for [/ ;^ t, A, whence the latter de¬ 
creases with increase in U. This shows clearly that IPT 
does not capture the spin physics of extreme correlations 
correctly. 

Within CTQMC, as is clear from Fig. ITOlfor A = l.Ot, 
as U/t increases, first the transition temperature Taf in¬ 
creases with increase in U/t for U/t < 5. This trend is 
similar to what is seen within IPT. But as Ujt increases 
further, the physics of effective Heisenberg model starts 
playing a role and the transition temperature starts de¬ 
creasing with further increase in C7 as it is governed by J. 
For [/ ;#> 2A, as A increases, the spin-exchange coupling 
J increases which is reflected clearly in the behaviour of 
Tat in Fig. [TTJ These results are consistent with earlier 
DMFT-I-CTQMC work [1^. Fig.fTTlshows the transition 
temperature Tjv as a function of A/t for a few values of 
U/t. We have shown comparison of T^r obtained within 
CTQMC with that of the Heisenberg model with spin 
exchange coupling of J. For U ^ 2A, J/4 is a very good 
approximation to Tat. But for U ^ 2A, H^ff is not the 



A/t 

FIG. 11: Top: Neel temperature Tn vs A/t obtained within 
DMFT-I-CTQMC for the IHM at half filling for various values 
of U/t. We see that for A -C C, Tn follows J/4 very closely. 
But for U ~ 2A, Tn starts deviating from J and decreases 
with increase in A/C/. 


correct low energy Hamiltonian of the model and we do 
not expect Tjv to be given by J/4. In fact in Fig. [TTl 
we see that Tm decreases as A increases beyond C//2 in 
contrast to what one would get from J/4. 

Finally, we present the full magnetic phase diagram 
within DMFT-I-CTQMC inferred from calculations for a 
range of parameter values in Fig. [121 As shown in the 
bottom panel of Fig[T3 in the 3 dimensional T — U — A 
space, there is a surface of first order phase transitions 
from PM BI to AFM insulator. Also there is a surface of 
second order phase transition across which the AFM or¬ 
der is lost continuously (although, as we have noted, this 
surface does not show up in IPT). These two surfaces 
are separated by a line of tri-critical points. This can be 
seen more clearly in the top panel of Fig. [121 Here the 
left panel shows mg vs U/t for various values of T. As 
T increases, the value of Uaf corresponding to the first 
order transition, where the AFM turns on with a jump, 
increases. This is because there are more thermal fluc¬ 
tuations and a larger U is required to stabilize the AFM 
order. Also, for the same reason, the AFM order does 
not survive for very small values of J and thus the Un 
at which the AFM order is lost by a continuous transi¬ 
tion decreases. These two transition points, namely Uaf 
(point of first order phase transition) and U at (point of 
second order phase transition) come close as T increases. 
There is a tricritical point which separates the two lines 
of first order and second order transitions. For A = l.Ot, 
from the CTQMC data we have generated, the tri-critical 
point seems to lie on the top of the dome of AFM region 
shown by a black point in the top-right panel of Fig. [T3 
but to be certain about this the calculations need to be 
done on a finer mesh of U/t values. 
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FIG. 12: Top: Left panel shows ms vs Ujt for A = l.Of 
and various values of /3. These results are obtained within 
DMFT+CTQMC. The mg turns on via a first order transition 
at Uaf{T) (shown as blue points in the right panel) while 
it is lost continuously at t/jv(r) (shown as red points in the 
right panel). As T increases, the U range, Un{T) — Uaf{T), 
in which the system shows AFM order shrinks to zero. A 
tri-critical point, shown as a black point in the right panel, 
separates the lines of first and second order phase transitions. 
Bottom: Phase diagram for the IHM at half filling in T — 
U — A space obtained within CTQMC. The surface made by 
the points connected by full lines is the first order transition 
surface from the PM to the AFM phase and the surface made 
by the points connected by dash lines is the second order 
transition surface from the AFM to a PM phase. 


DISCUSSION AND CONCLUSION 

In conclusion, in this paper we have presented sev¬ 
eral new results from a DMFT study of the ionic Hub¬ 
bard model at half filling, i.e., the Hubbard model in 
the presence of a staggered potential, which makes the 
system a HI for U = 0. As we turn on the on site re¬ 
pulsion U in this HI, first an AFM order sets in via a 
first order transition at 17 = Uaf- This is followed by 
a quantum phase transition to novel half-metallic AFM 
phase at 1/ = Uhm > Uaf- For still larger values of U, 
this system becomes an AFM insulator. Up to moder¬ 
ately strong values of U (e.g., U/t = 6.0 for A = l.Ot), 
the IPT captures the effects of electron-electron correla¬ 
tions quite well, and yields essentially the same results as 
CTQMC. But in the extremely correlated regime, where 
U ^ A, t, DMFT-blPT does not work well, as becomes 
clear when one does a finite temperature study. At any 
finite T, while the IPT continues to show only one first 


order phase transition at which the AFM order turns on, 
the CTQMC shows, in addition, a second, continuous 
transition back to a PM phase, with its physics deter¬ 
mined by the Heisenberg model. As T increases, the val¬ 
ues of U corresponding to the first and the second order 
transitions approach each other, shrinking the U range 
for which the long range AFM order is stable. There is 
a line of tricritical point Ttcp that separates the two sur¬ 
faces of first and second order phase transitions. To the 
best of our knowledge, this feature of the IHM has not 
been discussed earlier in the literature. 

Recently there has been a DMFT-bCTQMC study 
of the half hlled IHM within the PM sector, which shows 
a hrst order phase transition between Mott-Insulator and 
Metallic phases terminating at a critical point, just as in 
the Hubbard model at half Filing. However, this critical 
point lies inside the dome of the AFM region shown in 
top right panel of Fig. [12] and will be realised only if the 
AFM order is suppressed, either by lowered dimension¬ 
ality (eg., quasi 2-d systems) or due to frustration (eg., 
by the presence of next nearest neighbour hopping, or 
a frustrated lattice). We hope to study these issues in 
future work. At the end we would like to mention that 
recently the IHM has been realised in ultracold fermions 
on a 2-dimensional honeycomb lattice and it can be 
extended to higher dimensional layered honeycomb lat¬ 
tice by introducing perpendicular hopping. Though our 
numerical study is on the Bethe lattice of inhnite connec¬ 
tivity, we expect the qualitative physics to be the same 
for any bipartite lattice in d > 2 which has a compact 
density of states (DOS) like the DOS of the Bethe lattice 
of infinite connectivity. By choosing a large enough A, it 
might be possible to realise an AFM phase for the IHM in 
experiments where the AFM order turns on with a first 
order transition and is lost by a second order transition 
by tuning U. It would be interesting to look for signa¬ 
tures of the various effects we have discussed, including 
the quantum phase transition, in the experimental mea¬ 
surements in such systems. 
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U/t 


FIG. 13: Staggered magnetization rus plotted as a function of 
U/t for A = l.Ot. As pointed by the arrows, the red curve is 
obtained by doing the DMFT+IPT calculation for increasing 
U/t while the blue curve is obtained by starting from large 
U/t side and decreasing U/t. 


APPENDIX A 

To characterize the nature of a phase transition, one 
normally studies its hysteresis behaviour. We have solved 
the DMFT+IPT self-consistent equations, first starting 
from small U value and increasing U (AF-I) and then 
starting from large U value and decreasing U (AF-II). 
We see a clear hysteresis in the behavior of staggered 
magnetization mg and Sn. In Fig. 1131 we have shown 
results for the staggered magnetization for A = l.Ot. We 
see that in the AF-I solution, mg becomes non-zero for 
U > 3.0t. On the other hand, in AF-II solution, rUg re¬ 
mains non-zero up to t/ = 3.2t. To get the transition 
point Uaf, we compare the ground state energy in the 
PM phase with that in the AF-I solution and AF-II solu¬ 
tion. Fig, im shows the ground state energy for A = l.Ot 
as a function of U/t. For U < 3.0t, the PM phase is 
stable. For U > 3.0t, Egnd for the AF-I sector becomes 
lower than the ground state energy in the PM phase. No¬ 
tice that the Egnd of the AF-II sector becomes lower than 
the Egnd of the PM sector for larger value of U/t. Thus 
the AFM state becomes stable when for the first time mg 
becomes non-zero coming from the small U side. We call 
this point Uaf, which gives the boundary between PM 
BI and AFM phase in Fig. [2] 


APPENDIX B 

For the model in Eq. [1] of the paper, the self energy 
within the Hartree-Fock (HF) approximation is given by 

= U {uaa) = ^[^-Sn + anig] 

^B,a = U{nBi) = ^[l + 6n-amg] ( 22 ) 



U/t 

FIG. 14: The ground state energy Egnd plotted as a function 
of U/t for A = l.Ot. The red curve is Egnd obtained by solving 
the DMFT+IPT equation for model in Eq. ((T| in the spin- 
symmetric sector while the blue and the green curves show 
the Egnd obtained by solving the DMFT+IPT equations in 
the AFM sector. The blue curve is obtained by solving the 
self consistent equations coming from the small U side and 
the green one is obtained in decreasing U order. The AFM 
order sets in for U > Uaf = 3.Of for which the Egnd of the 
AFM sector is lower than that for the PM sector. 


Here nig = {mzA — nizB)/‘2 is the staggered magnetisa¬ 
tion with niza = tic-f- — Uai and a = A, B is the sublattice 
index. Sn = {nB — nA)/‘i^ is the staggered occupancy, i.e., 
the difference in the filling factor of the two sublattices. 
Since the bare Green’s function (with U=0) is given by. 


Gocr(k,iw„) = 


- A- 

-Ck 


-Ck 

- A- 


(23) 


the HF corrected Green’s function is given by 


Ga (k, iuJn) 


iUJn + ga +u -Ck 

-Ck iUJn - ga+ U 


-1 


(24) 


Here = ■^=0is the chemical potential and g^ = 

A — ^{6n + anig) which gives a gap Eger = \gcr\ in the 
single particle spectrum of a spin component. Using this 
Green’s function, one gets the following self consistent 
equations for the physical quantities defined above: 


Ws = ^ y dtpo{t)^^^[f{E,r{f)) - /(-A„(e))](25) 

= depo{e)^-^^[f[-Err{e)) - f{Eer{e))](2U) 


Here Ea{e) = + gl, f{Ea{e)) = exp(/3(i;+D-A))+i 

the Fermi function and po(e) is the bare density of states 
of the lattice under consideration. 
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FIG. 15: Top: Plots of ms and 5n vs U/t for A = l.Ot and 
n = 1 within the HF theory. Bottom: Spectral gap Eg^ for 
the up and down spin components within the HF theory. As 
soon as the magnetic order turns on, both Eg^f and Egi start 
increasing with U/t. Thus the half-metal phase, seen in the 
DMFT calculation just after the onset of the AFM order, 
is missing here and the system is an AFM insulator for all 
U > Uaf. 


We have solved the self-consistent equations for the 
Bethe lattice of infinite connectivity and the results ob¬ 
tained at half-filling (n = 1) and zero temperature are as 
follows. For small U/t the system is a BI with rus = 0 
and a non zero Sn. At U = Uaf a first order phase 
transition takes place with a jump in mg to a non zero 
value as shown in Fig. [151 For U > Uaf, the system 
is an AFM insulator. Fig. [15] shows the phase diagram 
at half-filling within the HF theory. For comparison we 
have also shown the phase diagram within DMFT+IPT 
at half-filling. The threshold Uaf required to turn on the 
magnetisation is smaller in the HF theory as compared to 
its value within the DMFT-blPT. This is because quan¬ 
tum fluctuations captured in DMFT are missing in the 
HF theory; as an effect the magnetic order survives up 
to smaller values of U. However, the HF phase transition 
line approaches the DMFT line as A/t gets smaller. The 
bottom panel of Fig.lTKlshows the spectral gaps Eg^. Just 
after the AFM order sets in, the spectral gaps for both 
the spin components start increasing with U/t which is 
in contrast to what is seen in the DMFT calculation. 



A/t 


FIG. 16: Phase diagram at half filling at T = 0 for Bethe lat¬ 
tice of infinite connectivity. Red circles are the data obtained 
from DMFT-I-IPT study while the blue circles are the data 
obtained from the HF theory. 


where the gap for one of the spin component keeps de¬ 
creasing with U/t even for U > Uaf leading to a HM 
point at Uhm > Uaf- Thus within the simple HF the¬ 
ory, where the self energy is independent of w, there is 
no half-metallic phase at half filling. 


APPENDIX C 


The ground state energy within the HF theory is 


Egnd — Ei{ek) U {na^){nai,) 

k k (a=A,B 

(28) 

where, as before, E„[ek) = -b with = A — 
^^{5n-\- ams). The last term in Egnd can be re-expressed 
as ^{6n^ — ml). Following the Ginzburg-Landau (GL) 
theory, we do the Taylor series expansion of Egnd for 
small mg'. 


E, 


gnd 


Eq + mlE"{ms = 0) -b mlE"" [mg = 0) -b 


= Eo + -TO, 


-TO„ 


-TO, 


6 


..(29) 


Here E" is second derivative of Egnd and so on. To de¬ 
cide about the nature of the phase transition, it is suffi¬ 
cient to look at the signs of the coefficients a, b and c 221 . 
For 6, c > 0, if a > 0, to, = 0 is the only point of minima 
of the ground state energy. As a changes sign, the system 
undergoes a second order phase transition to the magnet¬ 
ically ordered phase with ml = ■^{—b+ — 4ac). For 

c > 0 and 6 < 0, we have a first order phase transition 
at & = —tt\Jca/i where the magnetisation mg changes 

/ Q \ 1/4 

discontinuously by the amount (^) 
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8 



8 


96(C//2)« 


6 


4+2 4^2 

cl{g)£{- — )-c2{g)JC{- — ) 
9 9 


(30) 

with cl{g) = 32 + S2g^ + — g^ and c2{g) = 16 + 

24(7^+ 5 ^ — g®. Here/C(a;) = — xsin'^(d)]~^/^d 0 is 

the complete elliptic integral of the first kind and £ (x) = 


6 is the complete elliptic integral of 
the second kind. Fig. fTTl shows the plots of GL coefficients 
a, b and c (obtained from Eq. [30|) vs g for a fixed value 
of U. As 5 —>■ 0, /C(——>■ 0 while £{—^) —>■ oo. 
Thus for g < 2t, which is the regime of interest, b is 
always negative. Thus following the GL approach [ 2 ^ 
we conclude that the transition from the PM B1 to the 
AFM phase in the half filled IHM is always of first order 
in nature, even for very small values of A/t. 


20000 
15000 
O10000 
5000 
0 

0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 17: Plots of the GL coefficients a, b and c vs g for U = 
2.Of. One can see that a changes sign as g increases while 
6 < 0 and c > 0 for all values of g. 



Expressions for the GL coefficients in the Taylor series 
expansion of the ground state energy in Eq. are given 
below, 

Eq = -2 Y.k E{^k) - 

f = c/ + 2 (f)^E.^P-i] 

! = 6(f)^E.p^[l-6 r^ + 5r4] 
f = 90(f E. jEikw [-1 + - 35E + 21r6] 

Here E{ek) = E^{ek)\ms=o = \/^i+ 9 '^ with g = 
= A - f and r = 

We have numerically calculated the coefficients a, b 
and c and found that for all values of A and U/t studied, 
c is always positive while b is always negative, a > 0 
ioi U < Ui and becomes negative for 17 > 17i where the 
value of Ui depends upon A/t. For the Bethe lattice 
of infinite connectivity, the integrals involved in the 
above equations can be done analytically and we get the 
following expressions for the GL coefficients: 


1 = U + 


iU/2f 


Ag£{- 


^)_4l^^/C(- 
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